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ABSTRACT 

We calculate the string tension, deconfinement transition temperature and bulk 
thermodynamic quantities of the SU{3) gauge theory using tree level and tadpole 
improved actions. Finite temperature calculations have been performed on lattices 
with temporal extent A^^ = 3 and 4. Compared to calculations with the standard 
Wilson action on this size lattices we observe a drastic reduction of the cut-off depen- 
dence of bulk thermodynamic observables at high temperatures. In order to test the 
influence of improvement on long-distance observables at we determine the ratio 
Tc/ \fa. For all actions, including the standard Wilson action, we find results which 
differ only little from each other. We do, however, observe an improved asymptotic 
scaling behaviour for the tadpole improved action compared to the Wilson and tree 
level improved actions. 



1 Introduction 



Tree level and tadpole improved actions have been shown to yield a substantial 
reduction of cut-off dependences in the calculation of thermodynamic properties of 
the SU{3) gauge theory [|], In the high (infinite) temperature limit this is quite 
evident already from a perturbative calculation of the pressure (p) or energy density 
(e). In this limit high momentum modes give the dominant contribution to these 
observables. An improvement of the discretization scheme for the Euclidean action 
at short distances will thus naturally lead to a better representation of the ideal gas, 
Stefan-Boltzmann law on lattices with finite temporal extent N^-. However, even at 
Tc the improved actions lead to a reduced cut-off dependence for some observables. 
Calculations of the surface tension or the latent heat of the first order deconfinement 
transition, for instance, show a strong reduction of the cut-off dependence on lattices 
with temporal extent Nr = 3 and 4. In this case it also has been found that the 
remaining cut-off dependence on coarse lattices is weaker for a tadpole improved 
action than for a tree- level improved action [0. 

At high temperature the variation of p/T^ or e/T^ with temperature is small. 
Both observables only slowly approach the ideal gas limit, which is consistent with 
the expectation that thermodynamics depends on a running coupling that varies 
logarithmically with temperature. The numerical simulations performed in a given 
discretization scheme at finite temperature are therefore particularly sensitive to a 
correct representation of the infinite temperature limit. An accurate determination 
of the temperature scale itself is, however, not needed for the observation of the 
dramatic improvement in approaching the continuum Stefan-Boltzmann limit at 
infinite temperature. This has been utilized in the calculations presented in [||] where 
the temperature scale has been fixed using an effective coupling scheme combined 
with the asymptotic form of the SU(3) /^-function 0. At temperatures close to 
the deconfinement transition, however, thermodynamic observables like p/T'^ vary 
rapidly. A comparison of the improvement achieved with different actions in this 
temperature regime requires an accurate determination of the temperature scale. For 
some fixed point actions this has been done by determining the critical temperature 
on lattices with different temporal extent . Using calculations of to set the scale 
has, however, the disadvantage that T/Tc gets to be known only at a few discrete 
values. We will use here the string tension to define a continuous temperature scale, 
Tj^fa^ for thermodynamic observables. This also is needed for the determination 
of thermodynamic observables like the energy density, which require the knowledge 
of the SUi^) /5-function, i.e. the derivative of the bare couphng with respect to the 
cut-off in the non-asymptotic regime. 

In this paper we analyze the heavy quark potential with tree level and tadpole 
improved actions. From the potential the string tension is extracted in order to 
define a temperature scale, T j^/o. This in turn is used in a study of the thermody- 
namics of the SU(3) gauge theory. In the next section we discuss our calculations 
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Table 1: Tadpole improvement factors for the action S**-^'^^ defined in Eq. ( p.ll) . 



of the heavy quark potential and the determination of the string tension. The ap- 
plication of the resulting /3-function for thermodynamic calculations is presented in 
Section 3. Section 4 contains our conclusions. 



2 String Tension 

In our calculations we use two different improved actions for the SU(3) gauge theory. 
These actions include in addition to the standard 1x1 Wilson loop an additional 
1 X 2 or 2 X 2 loop, respectively, 



S(l,2) 



-5^'''^ = y (-W''l{x)-^W^/,ix)) . (2.1) 



X,U>fl 



Here uq denotes a tadpole improvement factor which we have chosen to be defined 
through the plaquette expectation value 0, i.e. uq = (1 - {W^:l{x))Y^\ The tree 
level improved actions are obtained for mq = 1. The high temperature ideal gas 
limit for these actions has been analyzed previously and first calculations of 
the pressure using the action S**^^'^-* have been presented there. The analysis of the 
infinite temperature ideal gas limit suggests, in fact, that the (l,2)-action is superior 
to the (2,2)-action. Although the leading O(a^) corrections are eliminated in both 
cases, it turns out that the remaining higher order contributions are much smaller 
for the (l,2)-action. In order to calculate the string tension from the long distance 
part of the heavy quark potential we have performed simulations of the SU (3) gauge 
theory on lattices of size 16^ and 24^ at several values of the gauge coupling j3 = G/g"^. 
The heavy quark potential has been extracted from smeared Wilson loops following 
closely the smearing approach described in [0], i.e. Wilson loops are constructed 
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from smeared links which are obtained by iterating the replacement process 



U^{x) ^ U,{x) + 7 E Uu{x)U^{x + u)Ul{x + fi) (2.2) 

several times. Some tests have been performed to find optimal values for the smear- 
ing parameter 7 for the range of couplings explored here. The value 7 varies between 
0.2 and 0.8. This allows to achieve a reasonable large overlap with the ground state 
already after about 10 to 15 smearing steps|. The potential at distance R is then 
determined from the asymptotic behaviour of smeared Wilson loops, W{R,L), 

Typically, for each value of the gauge coupling we have analyzed Wilson loops on 
300 configurations which were separated by 10 updates performed with an over- 
relaxed heat bath algorithm (1 update = 4 over-relaxation steps followed by 1 heat 
bath step). In the case of the (l,2)-action we also investigated the effect of tadpole 
improvement. For this purpose we have first determined self-consistently the factor 
uq at a few values of the gauge coupling. These numbers are given in Table |l|. 
Spline interpolations of these values have then been used for simulations at other 
values of the gauge coupling. Calculations of Wilson loops have been performed 
at 10 to 15 values of the gauge coupling. We generally observe that the ratios of 
smeared Wilson loops become, within errors, independent of L already for rather 
small values, i.e. for L > Lmin — (2 — 5) for large R. The potential has then been 
obtained from a weighted average of the logarithm of ratios W{R, L)/W{R, L -|- 1) 
with L > Ljnin- These calculations have been performed on lattices of size 16^, 
except for the largest values of the gauge coupling where calculations have been 
performed on a 24"^ lattice. Similar results have been obtained with the tree level 
improved action S*^^'^-* on lattices of size 24^. The potentials have then been fitted 
to the ansatz 

V{R) = Vo + ^ + aR , (2.4) 

for distances R > -Rmin in order to extract the string tension. For large values of R, 
roughly Ra > l/4fm [§, the coefficient of the Coulomb-like term is expected to be 
determined by string fiuctuations, d = — 7r/12 Indeed, for P corresponding to 
a lattice spacing smaller than a ~ 0.1 3fm we find values for a which are consistent 
with this (Table |^). Since we are interested in the long distance behaviour of the 
potential we followed the strategy to fix the coefficient of the Coulomb-like term for 
the determination of the string tension. The minimal value of -Rmin was chosen to 
be the maximum of Ra = 0.25fm and R = 2\/2, so that distortions of rotational 
symmetry which are present at small distances are not large. This approach was 
compared to a systematic increase of Rmin until stable results have been obtained for 
a. Typically this was the case for R^m^ 1 / (2"\/^) i.e. R^i^a> 0.25fm. We then have 

^Details on the optimization of our smearing procedure can be found in Ref. |^]. 
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Tree- level improved (l,2)-action 






-Ro.25fm 


min 
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16^ 


4.300 


1.98 


2.828 


-0.318 (82) 


16^ 


4.400 


2.34 


2.828 


-0.335 (68) 


16^ 


4.600 


2.97 


3.0 


-0.290 (36) 


16^ 


4.800 


4.11 


4.0 


-0.229 (76) 




5.000 


4.79 


5.0 


-0.284 (30) 


Tadpole improved (l,2)-action 




/5 


-Ro.25fm 


/?* 

min 


a 


16^ 


4.600 


1.96 


2.828 


-0.310 (61) 


16^ 


4.800 


2.53 


2.828 


-0.262 (42) 


24^ 


5.100 


3.74 


4.0 


-0.288 (40) 


24^ 


5.250 


4.13 


4.24 


-0.316 (112) 


24^ 


5.400 


4.76 


5.0 


-0.212 (54) 


Tree level improved (2,2)-action 






-Ro.25fm 


p* 

min 


a 


24^ 


4.600 


1.98 


2.828 


-0.299 (45) 


24^ 


4.800 


2.76 


2.828 


-0.296 (21) 


24^ 


5.000 


3.70 


4.0 


-0.274 (22) 


24^ 


5.200 


4.85 


5.0 


-0.271 (47) 


24^ 


5.400 


6.43 


5.0 


-0.289 (29) 


24^ 


5.600 


8.10 


5.0 


-0.266 (22) 



Table 2: Coefficient of the Coulomb-like term obtained from 3 parameter ffis to 
the heavy quark potential which have been calculated with tree level and tadpole 
improved (l,2)-actions and tree level improved (2,2) action. -Ro.25fm is the distance 
in lattice units corresponding to a physical value of 0.25 fm and -R^in minimal 
distance used in the ffis to the potential for the determination of a. The values of 
alpha and the errors have been calculated in the same way as for the string tension 
according to the formula ( |2.5| ) 
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Figure 1: The difference between the calculated potentials and the best fit obtained 
with fixed a in units of \fa versus Shown are results for the tree level improved 

(2,2)-action (a) and the tadpole improved (l,2)-action (b) for different values of the 
gauge coupling. Tree level results for the latter action look similar. 
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further increased -Rmin and have averaged the resuhs for d obtained from several fits 
with -Rmin > -Rmin Order to minimize remaining distortion effects resuhing from 
missing rotational invariance of the potential at these distances. Since the results 
for a at distances larger than -Rmin agree within errors we calculated the mean 
value according to the fit formula of a constant, but modified the error formula by 
introducing a factor ^J~N which takes into account that the different values of a are 
strongly correlated. 



cr 



y-v 0"(-Rmin) j [ ^ ^ 



Aa 



mm — mm / \ iiini — mm 

1 

, -'V 

\ ^?miii>-fimin 



2 



y — ' — ■ (2.5) 

iV . ^. ^ (Aa(/2min))^ ' 



The resulting values for the string tension for the different actions are given in 
Table ^ and in Table |[ Also given there are the minimal values Lmin and -Rmin used 
to extract and fit the potential for a given value of (3. In the last column we show the 
values of a obtained from the 3 parameter fit of the potential (Eq. p.4|) . Except for 
(3 = 4.400 for the tree level improved (1,2) and f3 = 4.500 for the tree level improved 
(2,2) action all values of the string tensions obtained from 2 parameter fits agree 
within errors with the ones from 3 parameter fits. The of the fits is comparable. 

In Figure we show the difference between the calculated potential and the 
best fit obtained with a = —n/12. Results are shown for the tree level improved 
(2,2)-action and the tadpole improved (l,2)-action. All fits have been performed 
for distances R^/&> 0.5. We note that the agreement between the numerical data 
and the fit generally remains good at distances smaller than the fitting interval. For 
smaller values of P we observe, however, a scattering of the data for Ry/d' < 0.5 due 
to the lack of rotational invariance. This effect becomes smaller for larger values of 
(3, i.e. closer to the continuum limit. For large values of (3 the potential can also be 
obtained at shorter physical distances. Here we clearly observe that the potential 
at small distances would prefer a larger coupling for the Coulomb term than used 
in our fit. This is reflected in the increase of Vdata — Vat at short distances. 



2.1 Deviations from asymptotic scaling 

Improved actions depend on several couplings for the different Wilson loops ap- 
pearing in the action. Through the specific choice of these couplings a particular 
trajectory in a mult i- dimensional parameter space is defined on which the continuum 
limit is approached. Although the improvement scheme for the actions does not aim 
at improving the approach to asymptotic scaling we may test in how far the scaling 
behaviour is modified through the specific choice of a trajectory. The results for the 
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Tree-level improved (l,2)-action 
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U. i ( O^Z iU^ J 


u. i ( Zoo ^^^o ( J 


16^ 


4.100 


3 


2.828 


0.14234 (165) 


0.13871 (842) 


16^ 


4.150 


3 


2.828 


0.11423 (131) 


0.11079 (586) 


16^ 


4.300 


4 


2.828 


0.07265 (107) 


0.06987 (489) 


16^ 


4.400 


4 


2.828 


0.05412 (51) 


0.05055 (182) 


16^ 


4.600 


4 


3.0 


0.03201 (43) 


0.03105 (182) 


16^ 


4.800 


5 


4.0 


0.01634 (48) 


0.01734 (257) 




5.000 


6 


5.0 


0.01188 (18) 


0.01104 (84) 


Tadpole improved (l,2)-action 




f3 


^min 


R* 

min 


0'2par 


''"Spar 


16^ 


4.050 


1 


2.828 


0.42159 (217) 


0.42216 (1462) 


16^ 


4.150 


2 


2.828 


0.29927 (340) 


0.32253 (3625) 


16^ 


4.185 


2 


2.828 


0.27153 (216) 


0.25530 (1645) 


16^ 


4.200 


2 


2.828 


0.26205 (198) 


0.25367 (1124) 


16^ 


4.250 


2 


2.828 


0.21829 (176) 


0.21355 (893) 


16^ 


4.300 


2 


2.828 


0.18343 (103) 


0.17852 (496) 


16^ 


4.350 


3 


2.828 


0.15416 (238) 


0.15947 (1265) 


16^ 


4.400 


3 


2.828 


0.12994 (154) 


0.12581 (754) 


16^ 


4.500 


3 


2.828 


0.09633 (53) 


0.09536 (252) 


16^ 


4.600 


4 


2.828 


0.07154 (90) 


0.06839 (366) 


16^ 


4.800 


4 


2.828 


0.04315 (55) 


0.04296 (235) 


24^ 


5.100 


6 


4.0 


0.01967 (22) 


0.01906 (87) 


24^ 


5.250 


6 


4.24 


0.01604 (53) 


0.01513 (233) 


24^ 


5.400 


6 


5.0 


0.01227 (20) 


0.01308 (97) 



Table 3: String tension obtained from 2 and 3 parameter fits to the heavy quark 
potential which have been calculated with tree level and tadpole improved (1,2)- 
actions. Lmin is the minimal extent of smeared Wilson loop W{R, L) used to extract 
the potential at distance R. -RJ^in ^^e minimal distance used in fits to the potential 
for the determination of the string tension. 
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Tree-level improved (2,2)-action 
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D* 

-'''min 






4.300 


3 


2.828 


0.24547 (981) 


0.31077 (7679) 


4.380 


3 


2.828 


0.17401 (278) 


0.15589 (1534) 


4.396 


3 


2.828 


0.16334 (228) 


0.16059 (1366) 


4.400 


3 


2.828 


0.16052 (207) 


0.16778 (1474) 


4.406 


3 


2.828 


0.15603 (186) 


0.14942 (1257) 


4.430 


3 


2.828 


0.13882 (141) 


0.13238 (900) 


4.450 


4 


2.828 


0.12787 (102) 


0.12698 (1949) 


4.500 


4 


2.828 


0.10422 (53) 


0.08755 (867) 


4.550 


4 


2.828 


0.08636 (83) 


0.08431 (461) 


4.600 


4 


2.828 


0.07060 (48) 


0.06957 (236) 


4.800 


4 


2.828 


0.03606 (13) 


0.03536 (57) 


5.000 


4 


4.0 


0.02046 (13) 


0.02021 (48) 


5.200 


5 


5.0 


0.01169 (16) 


0.01166 (76) 


5.400 


5 


7.0(5.0) 


0.00669 (14) 


0.00642 (46) 


5.600 


5 


8.0(5.0) 


0.00426 (14) 


0.00425 (35) 



Table 4: String tension obtained from 2 and 3 parameter fits to the heavy quark 
potential which have been calculated with the tree level improved (2,2)-action on 
lattices of size 24'*. L^in and -Rmin ^-re explained in Table |^. For the last two /3 
values we used a smaller value of -Rmin — 5-00 in the 3 parameter fit to obtain stable 
results. 
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Figure 2: y/a/ versus y/aa for various actions. Shown are results for y/aa/ R{g^) , 
normalized to Aj^. Here -R(5'^) denotes the 2-loop /5-function. In the case of the 
(2,2)-action the ratio of A^reeV^Ms known to us and we have, therefore, used 

an arbitrary normalization factor of 0.15. The horizontal line indicates the value of 
the continuum extrapolation of y/a/Ajj^ for the Wilson action taken from Ref. [6]. 



string tension given in Table |^ and | can be used to analyze the relation between the 
bare gauge coupling and the lattice cut-off, (3 (a) = 6/g'^{a). In the case of the Wilson 
action this has been analyzed in quite some detail [|T^. In a most straightforward 
way the deviations from asymptotic scaling become visible when one divides the 
results obtained for -y/aa by the universal 2-loop form of the /3-function. This yields 
^/a/AL. The perturbative expansion of the tadpole improved (l,2)-action coincides 
up to 0{g'^) with that of the tree level improved (l,2)-action with a modified gauge 
coupling /3 = /?(1 + a2C/V3), where = -0.366263(A^2 _ ly/^^ denotes the 0{g'^) 
expansion coefficient for the plaquette expectation value calculated with the tree 
level improved action |]l|. From this we obtain the ratio of A-parameters for tree 
level and tadpole improved (l,2)-actions, 

a(1,2) 

tad _ a2/6feo (n a\ 

A (1,2)-^ ' ^^-^l 



with 6o = llA^/487r . Using also the ratios of lattice A-parameters to ||Tl] , [12 
we obtain ^Jo j h.y^ for different actions. This is shown in Figure H. The numbers are 
plotted versus \foa in order to be able to compare results at the same value of the 
cut-off. The slope of these curves is an indication for the deviations from asymptotic 
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scaling in the regime of couplings investigated. We note that the tree level improved 
(2,2)-action shows similar scaling violations as the Wilson action while these are 
reduced for the tree level and even more for the tadpole improved (l,2)-actions. A 
similar behaviour is also obtained from calculations of the critical temperature with 
the tree level improved (1,2) action [|13 . 



3 Thermodynamics 

3.1 Tel \fo 

The string tension calculated in the previous section can be used to fix the tem- 
perature scale T I \fo = 1/ {Nr\fd). In particular, we have calculated the ratio 
Tc/ by determining the critical couplings for the different tree level and tadpole 
improved actions on lattices with temporal extent A^t = 3 and 4. Pseudo-critical 
couplings on finite spatial lattices have been determined from the location of the 
peak in the Polyakov loop susceptibility f^. For the tree level and tadpole improved 
(l,2)-actions we, furthermore, have performed a detailed study of the finite volume 
dependence of the critical couplings on lattices with temporal extent A^^ = 4 and 
N„ = 16, 24 and 32 0. In these cases the critical couplings have been extrapolated 
to the infinite volume limit using the ansatz 

/3,(Ar,, N^) = (3,iNr, cx)) - hi^-^j (3.1) 

which is appropriate for first order phase transitions. We note that in the continuum 
limit the difference Pc{Nt-,oo) — /5c(A'r, A'o-) is proportional to the finite volume 
shift in the critical temperature, i.e. (T^^oo ~ ^c,v)/^c,oo- In this limit Eq. ( |3.lD 
thus characterizes a physical property of QCD in a finite volume, the finite size 
dependence of peaks in susceptibilities which develop into a singularity in the infinite 
volume limit. Eq. ( p.l| ) then reads 

Tc,oo - Tcy ^ 47r2 h 

The parameter h thus is independent of the lattice action up to finite cut-off cor- 
rections. Given the current uncertainties in the numerical values of h these finite 
cut-off corrections can, however, not be disentangled from the statistical errors. 
We find for the tree level and tadpole improved (l,2)-actions proportionality factors 
h which are consistent with each other as well as with earlier results obtained with 
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other actions, 



h 



' 0.082 (32) 
0.072 (77) 
0.101 (34) 
0.068 (45) 
0.122 (54) 

^0.133 (63) 



;i,l) Wilson action [15], A^"^ = 4 
;i,l) Wilson action p[, A^^ = 6 
^1,2) tree level action, A^^ = 4 
^1,2) tadpole action, A^,- = 4 



(3.3) 



RG action 
RG action 



m 



3 
4 



The calculated critical couplings as well as the extrapolations to the infinite volume 
limit are summarized in Table § for the different actions. We note that our results 
for the tree level improved (l,2)-action are consistent with those of Ref. |]T3|. In 
Table ^ we also give results from [13| for larger values of A^,- and use the ansatz 
of Eq. ( ^.11 ) to extrapolate to the critical couplings on an infinite lattice. For the 
constant h we use in these cases the weighted average of the values given in Eq. ( p.3|) , 
i.e. h = 0.093 (18). 

In order to extract the critical temperature in units of the square root of the string 
tension we determine a at Pd^r, oo) from an interpolation with an exponential 
ansatz, a = Aexp{—P/2Nbo + f{P)), where f{j3) is a third order polynomial in 
The results for Tc/ \fa are shown in Table ^ and in Figure ^. 

The Wilson action /5c(A^r,oo) for =4 and 6 are taken from |15|. For A^^ = 8 and 
A^T- = 12 we performed the infinite volume extrapolation based on Eq. ( |3.1| ) with 
h = 0.093 by using j3c from 0. Recently there has been a new measurement of the 
string tension for the Wilson gauge action by Edwards et al. [0]. We used their 
"best" parameterization of the string tension to determine \/& at /3c- The results 
are displayed in Table |] and Figure |[ The strong cut-off effect so far reported for 
Tc/ ^/o^ is no longer visible in these new results For the continuum extrapolation 
we find the following value 



0.630 ±0.005 



(3.4) 



which is consistent with all ratios Tc/ ^/o^ extracted by us with improved actions on 
A^T- > 4 lattices (see Table |^). In Figure |^ and Table |^ we give the results obtained 
from calculations with different actions. This also includes results obtained with a 
renormalization group improved action [Q. The latter does lead to a slightly larger 
value for Tc/ ^/o^. However, also in this case the cut-off dependence is comparable to 
the other actions. This suggests that the difference between the results obtained with 
the RG-improved action in Ref. |TB[ and the results presented here is mainly due to 
differences in the analysis of the heavy quark potential rather than due to differences 
in the improvement scheme. In general we find that the cut-off dependence in the 



'^Previous estimates of Tc/ \/a- for the Wilson action at low p were based on preliminary results 
for the string tension |^ which propagated through the literature. Published values and 
recent high statistics results for the string tension are lower than the preliminary ones and 
thus lead to the increase in Tel \fa at A^,- — 4 and 6. 
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Tree level improved (l,2)-action 


a 






Tc/ a/o" 




3 


3.9079 (6) 




{oof 


3 


3.9094 (6)(3) 


0.630 (5) 


32^ 


4 


4.0729 (3) 




{oof 


4 


4.0730 (3) 


0.636 (4) 


20^ 


5 


4.19963 (14) 




{oof 


5 


4.20108 (14)(28) 


0.631 (5) 


24^ 


6 


4.31466 (24) 




{oof 


6 


4.31611 (24)(28) 


0.632 (5) 


Tree level improved (2,2)-action 










24^ 


4 


4.3995 (2) 




{oof 


4 


4.3999 (2)(1) 


0.625 (4) 


Tadpole improved (l,2)-action 






/3c 




12^ 


3 


4.1868 (4) 




{oof 


3 


4.1882 (4)(3) 


0.643 (3) 


32^ 


4 


4.3522 (4) 




{oof 


4 


4.3523 (4) 


0.639 (6) 



Table 5: Critical couplings for the tree level and tadpole improved actions. In each 
case we give the result for the largest spatial lattice on which simulations have been 
performed and the infinite volume extrapolation. Details on the determination of 
critical couplings for the (l,2)-actions for A^,- = 3 and 4 on lattices with different 
spatial extent are given in [2]. The finite lattice results for A^^- = 5 and 6 are taken 
from [13]. Infinite volume extrapolations are based on Eq. (|3.1| ) with h = 0.093, 
except for the case of the = 4 (l,2)-actions where a detailed finite volume scaling 
analysis has been performed in [2]. The second error on jSc in the infinite volume 
limit is the systematic error due to the error on h. In the last column we also give 
the ratio Tc/ \fa. The errors on Tc/ i/a due to the error on /5c and on \fd have been 
added. 
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Wilson action 


ATS 




/3c 




(oo)3 


4 


5.6925 (2) 


0.627 (6) 


(cx))3 


6 


5.8941 (5) 


0.632 (11) 


(oo)3 


8 


6.0624 (9) (3) 


0.629 (6) 


(cx))3 


12 


6.3380 (13) (10) 


0.630 (5) 



Table 6: Critical couplings for the Wilson gauge action. We give the results for 
the infinite volume extrapolations for A^^ = 4 and 6 from [15] and for A^^ = 8 and 
12 based on Eq. ( p.l|) with h = 0.093, where the values for /5c on finite lattices are 
taken from [3] . The second error given is the systematic error due to the uncertainty 
in the extrapolation. In the last column we give the ratio Tc/ ^/a using the string 
tension parameterization given in [14]. 



0.68 



0.66 



0.64 




0.62 



0.58 



Figure 3: The critical temperature in units of the square root of the string tension 
for various actions versus the square of the cut-off. The iV^ = 6 point of the tree 
level (1,2) improved action has been slightly shifted to make it distinguishable from 
the Wilson action point. 
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ratio Tc/ ^Jo is quite small for all actions. We also note that the results on A^,- = 3 
lattices are in good agreement with those obtained on larger lattices. This is quite 
different from the surface tension analysis It may, however, indicate that the 
latter is more sensitive to high momentum modes than the ratio Tc/a/^, which is 
also reflected in the fact that the Wilson action does not show large finite cut-off 
effects in T^/ s/a whereas there are large finite cut-off dependences in the surface 
tension. 



action 


/9c 




standard Wilson 
(2,2) (tree level improved) 
(1,2) (tree level improved) 
(1,2) (tadpole improved) 
(1,2) (RG improved) 


5.69254 (24) 
4.3999 (3) 
4.0730 (3) 
4.3523 (4) 
2.2879 (11) 


0.627 (6) 
0.625 (4) 
0.636 (4) 
0.639 (6) 
0.653 (6)(1) 



Table 7: Critical temperature in units of ^ on lattices with temporal extent A^,- = 4. 
Infinite volume extrapolations for the critical couplings have been performed in all 
cases. Further details on the data for the RG-improved action can be found in [18]. 



3.2 Pressure of the SU(3) gauge theory 

The temperature dependence of the pressure can be obtained from an integration 
of action densities, (S"), calculated at zero and non-zero temperature, respectively^. 



2^4 



' = K\yfi\{S\-{S)T) . (3.5) 



Here /3o is a value of the coupling constant below the phase transition point at 
which the pressure can safely be approximated by zero. The subscripts and T 
refer to calculations of the action expectation values on the zero temperature and 
non-zero temperature lattices, respectively. We also note that in the case of tadpole 
improved actions the action itself depends implicitly on the gauge coupling through 
the tadpole factor mo(/3). This has to be taken into account in the calculation of 
derivatives with respect to /3. We therefore have introduced in Eq. ( p.5|) the quantity 
S defined as 

S = S-/3|. (3.6) 

Using the string tension values given in Table ^ and ^ and normalizing these to the 
string tension at /3c, we obtain the temperature in units of T^. The temperature 



^We refer to Refs. & M for more details on the formalism. 
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at intermediate values of the coupling has been obtained from an interpolation. 
With this we can reanalyze also the results for the pressure obtained for the tree 
level improved (2,2)-action given in [|l|. It simply amounts to a modification of the 
temperature scale, i.e. the ordinate of Figure 3 in [^. While this has practically no 
consequences at high temperature, it leads to visible shifts close to T^. In Figure | 
we show the results of calculations with improved actions on 16^ x 4 lattices. These 
are compared to the continuum extrapolation obtained from calculations with the 
Wilson action on lattices with temporal extent A^^ = 4, 6 and 8 0. We find that 
all improved action calculations are close to the continuum extrapolation, while the 
standard Wilson calculation on a A^^ = 4 lattice clearly deviates substantially. 

We furthermore have calculated the pressure for the tree level and tadpole im- 
proved (l,2)-actions on a lattice of size 12^ x 3. These results are shown in Figure 
Here we also show the result of a calculation using a fixed point action for Nt = 3 
1^. These are in good agreement with the continuum extrapolation obtained from 
the Wilson action and seem to be closer to the result obtained with the tadpole 
improved action than the tree level result. For temperatures in the range (2 — 3)Tc 
the tadpole action yields about 10% smaller values for the pressure than the tree 
level action, although both approach the same infinite temperature limit. We also 
note that we do not observe any significant cut-off dependence when comparing cal- 
culations on lattices with temporal extent A^^ = 3 and 4 despite the fact that the 
cut-off dependence in the infinite temperature limit leads to about 15% differences 
in the Stefan-Boltzmann limit. This is, to some extent, in accordance with the 
analysis of the Wilson action, where we noted already that the cut-off distortion in 
this temperature range is only half as large as expected on the basis of calculations 
for the ideal gas limit It seems that these are further reduced in calculations 
with improved actions. Using the temperature scale defined by the string tension 
calculations we also can extract the /3-function, adj3/da, outside the validity regime 
of the asymptotic 2-loop form, i.e. in the coupling range explored here. With this, 
further thermodynamic observables can be extracted. In particular, we obtain for 
the difference between the energy density, e, and 3p, 



e-3p fKV'/ dp 
a 



T4 \N„J V da 



{S)o - {Sh 



(3.7) 



In Figure ^ we compare the result obtained with tree level and tadpole improved 
actions on lattices with temporal extent N^- = 4 with corresponding results obtained 
with the Wilson action for A^,- = 8. In the case of the Wilson action it has actually 
been observed that results for A^^ = 6 and 8 coincide within errors and may thus 
be taken as the continuum limit result [Q. The good agreement we find here with 
the Wilson action calculation confirms this observation, i.e. cut-off effects are indeed 
small in (e — 3p)/T^. In fact, the strong cut-off effects present in the ideal gas limit 
cancel exactly in this quantity, which in the high temperature limit is 0{g'^{T)). The 
good agreement between results obtained with different actions also is an excellent 
consistency check for the determination of the temperature scales and the analysis 
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Figure 4: Pressure of the SU(3) gauge theory calculated with the Wilson action 
and different improved actions on A^,- = 4 lattices (upper figure). The lower figure 
shows a comparison of calculations with tree level and tadpole improved actions on 
Nt- = 3 and 4 lattices. Also shown there are results from a calculation with a fixed 
point action (triangles). The arrows indicate the ideal gas result in the continuum 
limit. For comparison we also show the continuum extrapolation obtained from 
calculations with the Wilson action [3]. 
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Figure 5: The difference (e — 3p)/T^ versus T/Tc calculated on lattices with tem- 
poral extent A^^ = 4 for various improved actions. The results are compared to a 
calculation with the Wilson action on a A^^ = 8 lattice. In the lower figure the peak 
of (e — 3p)/T^ is shown in detail. 
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performed here. 



4 Conclusions 



We have analyzed thermodynamic properties of the SU(3) gauge theory using tree 
level and tadpole improved actions. In general we find that the use of improved 
actions does lead to a significant reduction of the cut-off dependence in observables 
like the energy density, the pressure and even the surface tension at T^. The im- 
provement has, however, little effect on the calculation of long-distance quantities 
like Tc/ On the other hand we do observe an improved scaling behaviour for the 
tadpole improved action. There also might be a slight advantage in the use of the 
tadpole improved action for the analysis of bulk thermodynamic observables like the 
pressure. However, at present this cannot be further quantified within the accuracy 
of our calculations. 

It now seems that the systematic cut-off dependencies in calculations of thermo- 
dynamic observables (T^, equation of state, latent heat, surface tension,..) are well 
under control with presently used improved actions. Remaining ultra-violet cut-off 
dependences, which without doubt still are present, are hidden by statistical errors 
and/or inaccuracies in the determination of the observables. The latter result, for 
instance, also from infra-red cut-off effects. This is true also for the determination 
of the string tension, which is sensitive to the specific form of the fit used to ana- 
lyze the heavy quark potential at distance R ~ (0.25 — 1) fm. Such ambiguities are 
likely to be the origin for the currently existing discrepancy between the calculations 



presented here and in |]T8|. However, also the thermodynamic calculations are still 
sensitive to infra-red effects. While the finite volume effects are quite well under 
control for the determination of the critical temperature they certainly have to be 
analyzed in more detail for the discontinuities (surface tension, latent heat) at T^. 
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